Spectra of complex networks 
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We propose a general approach to the description of spectra of complex networks. For the 
spectra of networks with uncorrelated vertices (and a local tree-like structure), exact equations 
are derived. These equations are generalized to the case of networks with correlations between 
neighboring vertices. The tail of the density of eigenvalues p (A) at large |A| is related to the behavior 
of the vertex degree distribution P (k) at large k. In particular, as P (k) ~ fc~ 7 , p (A) ~ |A| 1-27 . We 
propose a simple approximation, which enables us to calculate spectra of various graphs analytically. 
We analyse spectra of various complex networks and discuss the role of vertices of low degree. We 
show that spectra of locally tree-like random graphs may serve as a starting point in the analysis of 
spectral properties of real-world networks, e.g., of the Internet. 
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I. INTRODUCTION 

Many real-world technological, social and biological 
complex systems have a network structure. Due to their 
importance and influence on our life (recall, e.g., the In- 
ternet, the WWW, and genetic networks) investigations 
of properties of complex networks are attracting much 
attention [1-7]. Such properties as robustness against 
random damages and absence of the epidemic thresh- 
old in the so called "scale-free" networks are nontriv- 
ial consequences of their topological structure. Despite 
undoubted advances in uncovering the main important 
mechanisms, shaping the topology of complex networks, 
we are still far from complete understanding of all pecu- 
liarities of their topological structure. That is why it is 
so important to look for new approaches which can help 
us to reveal this structure. 

The structure of networks may be completely described 
by the associated adjacency matrices. The adjacency ma- 
trices of undirected graphs are symmetric matrices with 
matrix elements, equal to number of edges between the 
given vertices. The eigenvalues of an adjacency matrix 
are related to many basic topological invariants of net- 
works such as, for example, the diameter of a network 
[8,9]. Recently, in order to characterize networks, it was 
proposed to study spectra of eigenvalues of the adjacency 
matrices as a fingerprint of the networks [10-17]. The 
rich information about the topological structure and dif- 
fusion processes can be extracted from the spectral anal- 
ysis of the networks. Studies of spectral properties of 
the complex networks may also have a general theoret- 
ical interest. The random matrix theory has been suc- 
cessfully used to model statistical properties of complex 
classical and quantum systems such as complex nucleus, 
disordered conductors, chaotic quantum systems (see, for 
example, reviews [18]), the glassy relaxation [19] and so 
on. 

As the adjacency matrices are random, in the limit 



N — > oo (N is the total number of vertices), the den- 
sity of eigenvalues could be expected to converge to the 
semicircular distribution in accordance with the Wigner 
theorem [20]. However, Rodgers and Bray have demon- 
strated that the density of eigenvalues of a sparse random 
matrix deviates from the Wigner semicircular distribu- 
tion and has a tail at large eigenvalues [21], see also [22]. 
Recent numerical calculations of the spectral properties 
of small- world and scale-free networks [12-14], and the 
spectral analyses of the Internet [10,11,15,16] have also 
revealed that the Wigner theorem does not hold. The 
spectra of the Internet [10,11] and scale-free networks 
[13,14] demonstrate an unusual power-law tail in the re- 
gion of large eigenvalues. At the present time there is 
a fundamental lack in understanding of these anomalies. 
In order to carry out a complete spectral analysis of real 
networks it is necessary to take into account all features 
of these complex systems described by a degree distribu- 
tion, degree correlations, the statistics of loops, etc. At 
this time there is no regular approach that allows one to 
handle this problem. Our paper fills this gap. 

Our approach is valid for any network which has a 
local tree-like structure. In particular, these are uncor- 
related random graphs with a given degree distribution 
[23,24], and their straightforward generalizations [25] al- 
lowing pair correlations of the nearest neighbors. These 
graph ensembles have one common property: almost ev- 
ery finite connected subgraph of the infinite graph is a 
tree. The tree is a graph, which has no loops. A ran- 
dom Bethe lattice is an infinite random tree-like graph. 
All vertices on a Bathe lattice are statistically equivalent 
[26] . These features (the absence of loops and the statisti- 
cal equivalence of vertices) are decisive for our approach. 
The advantage of Bethe lattices is that they frequently 
allow analytical solutions for a number of problems: ran- 
dom walks, spectral problems, etc. 

Real- world networks, however, often contain numerous 
loops. In particular, this is reflected in a strong "cluster- 
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ing" , which means that the (relative) number of loops 
of length 3 do not vanish even in very large networks. 
Nevertheless, we believe, that the study of graphs with a 
local tree-like structure may serve as a starting point in 
the description of more complex network architectures. 

In the present paper we will derive exact equations 
which determine the spectra of infinite random uncorre- 
cted and correlated random tree-like graphs. For this, 
we use a method of random walks. We propose a method 
of an approximate solution of the equations. We shall 
show that the spectra of adjacency matrices of random 
tree- like graphs have a tail at large eigenvalues. In the 
case of a scale-free degree distribution, the density of 
eigenvalues has a power-law behavior. We will compare 
spectra of random tree-like graphs and spectra of real 
complex networks. The role of weakly connected vertices 
will also be discussed. 



II. GENERAL APPROACH 

Let A = (a vw ) be the N x N symmetric adjacency ma- 
trix of an ./V- vertex Mayer's graph G, a 2 uw = a vw , a vv = 
(the Mayer graph has either or 1 edges between any pair 
of vertices, and has no "tadpoles", i.e., edges attached at 
a single vertex) . Degree k v (the number of connections) 
of a vertex v is defined as 



k v — ^ ^ CLv 



(1) 



A random graph, which is, in fact, an ensemble of graphs, 
is characterized by a degree distribution P(k): 



I 1 N 



k) 



(2) 



Here, ( ) is the averaging over the ensemble. We suppose 
that each graph in the ensemble has N vertices. Graph 
ensembles with a given uncorrelated vertex degree dis- 
tribution may be realized, e.g., as follows. Consider all 
possible graphs with a sequence of the numbers {N (k)} 
of vertices of degree k, k — 1,2,..., J2k N(k) — N, as- 
suming N (k) /N — > P (k) in the thermodynamic limit 
[N(k) — > oo, TV — > oo]. Suppose that all these graphs are 
equiprobable. Then, simple statistical arguments lead to 
the conclusion that almost all finite connected subgraphs 
of an infinite graph do not contain loops. 

This approach can be easily generalized to net- 
works with correlations between nearest-neighbor ver- 
tices, characterized by the two- vertex degree distribution: 

P 2 (M') = (^- E a vw S(k v -k)S(k w -k')\ . (3) 



v,w — 1 



Here L = (1/2) J2 V w a vw is the total number of edges. 
In the case of an uncorrelated graph we have 



kk' 

P 2 (k,k') = —^P(k)P(k') 
(k) 



(4) 



where (k) = 2L/N is the mean degree of a vertex. 

The spectrum of A may be calculated by using the 
method of random walks on a tree-like graph G and gen- 
erating functions [27]. We define a generating function 
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(5) 



v—1 n—0 



where p v (n) is the number of walks of length n from v to 
v, where v is any vertex of G: 



p v {n) = {A n ) v , v . 



(6) 



In a tree-like graph the number of steps n is an even 
number. In order to return to v we must go back along 
all of the edges we have gone. 

Let q v (n) be the number of walks of length n starting 
at v and ending at v for the first time. We define 



Qv(z) = ^2q v {n)z n . 
One can prove that 
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(7) 



n=0 



N 



Qv(z) 



(8) 



Let d(w, v) — m ^ 1 be the distance from to to !i and 
iitvu (n) be the the number of paths of length n starting 
at w and ending at v for the first time. We define 



T^\z) = Y^t^ v {n)z n . 



n=0 



One can prove 

Q v (z)=zJ2T^(z), 

w 



(9) 



(10) 
(11) 



where w — ► gi — ► g 2 — ► •••<?m-i — > v is the shortest 
path from w to v. There is an important relationship: 



T<U(z) = z + zY i T§Kz) 
g 

= Z |z^T 9 W(z)tW(z). 



(12) 



In this sum the vertex g is the nearest neighbor of w and 
a second neighbor of the vertex v. Solving the recurrence 
equation (12), we can find TwJ(z) and Q v (z). 

We define T${z) = ^{z' 1 ). Equation (12) may be 
written in a form 



2 



(13) 



We can find Q v (z), from which we get R(z). Let us define 
B{z) = z^ 1 R{z^ 1 ). Then the density of the eigenvalues 
A of a random graph is determined as follows: 



p{\) = -Im(S(A + fe)))/7 



(14) 



where e is positive and tends to zero. Note that the 
equations (7)-(13) are valid for both uncorrelated and 
correlated tree-like graphs. 

In the case of a fc— regular connected graph we have 
f^(z) = T(z) and Q v (z) = Q{z). Eq. (13) gives 

zT(z)-(fc-l)T 2 (z) = l. (15) 

Solving this equation, we get the well known result: 



k \/4(fc — 1) - A 2 
2^ fc 2 -A 2 



(16) 



This is a continuous spectrum of extended eigenstates 
with eigenvalues |A| < 2y/k — 1. The presence of the de- 
nominator on the right-hand side of Eq. (16) leads to a 
difference of the spectrum of this graph from Wigncr's 
semi-circular law. In exact terms, Wigner's law is valid 
for the eigenvalue spectra of real symmetric random ma- 
trices whose elements are independent identically dis- 
tributed Gaussian variables [20]. These specific random 
matrices for Wigner's law essentially differ from the ad- 
jacency matrices, which we consider in this paper. So, 
in our case, the semicircular law may be used only as a 
landmark for a contrasting comparison. 



III. SPECTRA OF UNCORRELATED GRAPHS 

In the case of uncorrelated random tree-like graphs, 
k w — 1 random parameters Tgw(z) on the right-hand side 
of Eq. (13) are equivalent and statistically independent. 
They are also independent on the degree k w . We define 
the distribution function of T^v (z) at z = A + ie in the 
Fourier representation as: 



F x (x) = (exp 



cT^(X + ie) 



(17) 



where the brackets (...) means the averaging over the en- 
semble of random uncorrelated graphs associated with a 
degree distribution P(k). The statistical independence 
of the k — 1 random parameters Tgw(\ + ie) = Ti, 
i = 1,2, ...fc — 1, k = k w , on the right hand side of Eq. 
(13) allows us to use the following identity: 



F x {x) = (exp(-«T)) = ( exp(- 



IX 



) 



oo 

= i - v^J ^Ji(2v^y 



k-1 



(18) 



where J\{x) is the Bessel function and (fc) = J2k fcP(fc). 
Thus, we get the exact self-consistent equation for F\(x): 

F x (x) = 1 - y/E ^ ^Lj^^Eyyy^F^y)) , (19) 

Jo Vv 

where $i(x) = J^kLi kP^x^ 1 / (fc). Solving Eq. (19) 
gives the distribution of T, and so we can obtain Q, from 
which we get R. Eqs. (8), (10), and (14) give 



p(X) = Im 



* \A-E-=i^/ 
= -Rc/ dye^F^y)) , 

1" JO 



(20) 



where <f>(x) = T,T=i p ( k ) xk - From E( l- ( 19 )' wc nnd 
the n-th moment of the distribution function ^ X (T), Eq. 
(17): 

1 f°° 

M n ee (T n ) = {n _ 1)Hn J dyy^e^iF^y)) . 

(21) 



IV. EFFECTIVE MEDIUM APPROXIMATION 

In a general case it is difficult to solve Eq. (19) ex- 
actly. Let us find an approximate solution. We neglect 
fluctuations of T around a mean value T(A) = (T). A 
self-consistent equation for the function T(A) may be ob- 
tained if we insert 

F x (x) w e - txT ^ (22) 

into the right-hand side of Eq. (21) for n = 1. We get 



(fc) ^ A + ie - (fc - l)T(A) ■ 



(23) 



Below we will call this approach an "effective medium" 
(EM) approximation. At real A, T(A) is a complex func- 
tion, which is to be understood as an analytic continua- 
tion from the upper half-plane of A, T (A) = T (A + ie). 
Therefore, lmT(X + ie) < 0. In the framework of the EM 
approach, the density p(A), Eq. (20), takes an approxi- 
mate form 

fcP(fc) ImT(A) 



P(A) = ~E 



tt^ (A-fcReT(A)) 2 + fc 2 (ImT(A)) 2 



(24) 
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V. TAIL BEHAVIOR AND FINITE-SIZE EFFECT 

Equation (23) may be solved analytically at |A| ;» 
1. We look for a solution in the region ImT(A) <C 
RcT(A) « 1. It is convenient to use a continuum ap- 
proximation in Eq. (23). The real and imaginary parts 
of this equation take a form 



ReT(A) = 

rk cut 
Jko 

1 = 

/' 

/A "l-(fc. 



1 



2A (k) 



dkkP{k) 



l_(fc_l)Re^y 

1 



2 ' 



(25) 



A 2 (k) 



dkk{k-\)P{k) 



(26) 



where ko and k cut are the smallest and largest degrees, 
respectively. A region ko < k <C k\ gives a regular 
contribution into the integrals (25) and (26) while a re- 
gion k ~ k\ » 1 gives a singular contribution. Here 
k\ = A/ ReT(A) + 1. As a result we obtain 



2A 2(fc)ImT(A) ' 
A ' A -A/, A /';/••> 



l -^(k)l k X dkk{k - l)m + 



(k) ImT(A) 



(27) 
(28) 



If P(fc) decreases faster than fc 2 at k ^> 1, i.e. (fc) is 
finite, then in the leading order of 1/A wc find 



T(X) = \- 1 -iw\X\k x P(k x )/(k) . 



(29) 



Within the same approach one can find from Eq. (24) 
that the density p(X) also has two additive contributions 



(fc)ImT(A) k x P(k x ) 

Inserting Eq. (29) gives the density 
p(A)-2^P(fc A ). 



(30) 



(31) 



Here k x = A/ RcT(A) + 1 = A 2 + 0(1). 

The asymptotic expression (31) is our main result. The 
right-hand side of this expression originates from two 
equal, additive contributions: the contribution from the 
real part of T{\) and the one from the imaginary part of 
T(A). One can show that the asymptotic behavior of the 
real part, ReT(A) = A -1 + 0(A~ 3 ), in the leading order 
of 1/ A is universal and is valid even for graphs with finite 



loops. Contrastingly, the asymptotics of ImT(A) in the 
leading order of 1/ A and the corresponding contribution 
to the right-hand side of Eq. (31) depend on details of 
the structure of a network. 

The analysis of Eq. (23) shows that the main contribu- 
tion to an eigenstate with a large eigenvalue A is given by 
vertices with a large degree k ~ k\ » 1. As we shall show 
below, in the limit A » 1, the result (31) is asymptoti- 
cally exact. The relationship between largest eigenvalues 
and highest degrees, A 2 + O(l) = fc, for a wide class of 
graphs was obtained in a mathematical paper, Ref. [28]. 
This contribution of highly connected vertices may be 
compared with a simple spectrum of "stars", which are 
graphs consisted of a vertex of a degree k, connected to 
k dead ends. The spectrum consists of two eigenvalues 
A = ±^/k and a (k— l)-degenerate zero eigenvalue. Note 
that asymptotically, in the limit of large A, Eq. (31) gives 
p(X) = 2|A|P(A 2 ) if P(k) decreases slower than an expo- 
nent function at large k, that is, if higher moments of the 
degree distribution diverge. 

A classical random graph [29,30] has the Poisson de- 
gree distribution P(k) = e~ < fe > (k) k jk\. The tail of p(X) 
is given by Eq. (31) with k\ = A 2 + a 3> 1 where a is a 
number of the order of 1: 



p(X) ~ A- 2 < A2+a ) cxp[(l + ln(fc))A 2 ]. 



(32) 



This equation agrees with the previous results [21,22] ob- 
tained by different analytical methods. 

For a "scale-free" graph with P(k) ~ Pofc -7 at large 
k, at |A| S> 1, we get an asymptotically exact power- law 
behavior: 



p(A) £* 2 |A| P(A 2 ) = 2P |A| 



(33) 



where the eigenvalue exponent 5 = 2-y — 1. 

At a finite N ^> 1, there is a finite-size cutoff of the 
degree distribution k cut ~ /c ^V 1/(7_1) [31]. The cutoff 

1 /2 

determines the upper boundary of eigenvalues: A < k c ^ t . 
This result agrees with an estimation of the largest eigen- 
value of sparse random graphs obtained in Ref. [32]. 

Let us analyze the accuracy of the EM approach. One 
can use the following criterion. We introduce a quan- 
tity q n ee M„/T"(A). Here, T"(A) is the n-th mo- 
ment of the approximate distribution (22). Inserting 
the function (22) into Eq. (21) gives M n . The func- 
tion F\[x) = e~ lxT ( x ) would be an exact solution of Eq. 
(19) if q n = 1 for all n ^ 1. Note that at n = 1 wc 
have q\ = 1, because this equality is the basic equation 
in the framework of the EM approximation. At A > 1 
and P(k) = P /c~ 7 , in the leading order of 1/A, Eq. (21) 



gives 



q n = l- Ank \- 2 at 7 > 3 , (34) 

= 1 - Anfc A" 2 lnA at 7 = 3 , (35) 
= 1 - Ank^X- 2 ^-^ at 2 < 7 < 3 , (36) 
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where kg is the smallest degree in P{k) and A is a numer- 
ical factor. This estimation allows us to conclude that at 

1/2 

A S> k ' the EM solution becomes asymptotically exact. 

1 11 

At small A < k , the EM approximation are less ac- 
curate. For example, at A = for a scale-free network, 
we obtain q n = S n /S?, where S n = (k(k - l)~ n ) / (k). 
Only at large 7>1 and n < 7, the parameter q n is close 
to 1, i.e. q n = 1 + 0(n/7). 

One can conclude that the EM approach gives a reli- 
able result close to the exact one in the range 



fcy 2 «A«fc i ,/? = fcy 2 iv i /2(7-i). 



(37) 



In our derivations we assumed the tree-like local struc- 
ture of a network, that is, the absence of finite-size loops 
in an infinite network. Loosely speaking, this assumption 
may fail if the second moment (k 2 ) of the degree distribu- 
tion diverges. This can be seen from the following simple 
arguments. The length of a typical loop is of the order 
of the average shortest-path length of a network. Since 
the mean number of the second-nearest neighbors in the 
infinite uncorrelated net is (k(k — 1)) and diverges if (fc 2 ) 
diverges, the average shortest-path length and the length 
of a typical loop are small and may turn out to be finite 
even in the limit of an infinite net if (k 2 ) diverges. In 
this situation, the result (whether there are loops of fi- 
nite length in the infinite network or not) is determined 
by the size-dependence of the cut-off of the degree dis- 
tribution [24] . In its turn, this dependence is determined 
by the specifics of an ensemble and varies from network 
to network. 



VI. SPECTRA OF CORRELATED GRAPHS 

Many real-world networks are characterized by strong 
correlations between degrees of vertices [33-35] . The sim- 
plest ones are correlations between degrees of neighboring 
vertices. Let us study the effect of degree correlations on 
spectra of random tree-like graphs. 

Using the pair degree distribution (3), it is convenient 
to introduce the conditional distribution that a vertex of 
degree k is connected to a vertex of degree k\ : 



P(fci \k) = (k)P(k 1 ,k)/kP(k). 



(38) 



The method used above for the calculation of spectra 
of uncorrelated graphs may be generalized to correlated 
graphs. For this, one should take into account correla- 
tions between the degree of a vertex v and the generating 
function TyJJ(z) in Eq. (13). We define the distribution 
function of T wv in the Fourier representation as: 



F x (k, x) = (5 (k v - k) exp \-ixT^(\ + ie) 



(39) 



Averaging Eq. (13) and using the identity (18), we obtain 
an exact equation for F\(k,x): 



F x (k,x) = l-^^ l P(k 1 \k) 

fei 

x r^J 1 (2V^)e^ A F A fel - 1 (fc 1 , 2 ;), (40) 

Jo vy 

The density of eigenvalues is of the following form 

p(A) = -Re^P(fc) / dye^F^y). (41) 

These equations are a generalization of the equations 
derived above for uncorrelated graphs. Indeed, for an 
uncorrelated graph, we have P{k\ \k) = kP(k)/ (k) and 
F\(k,x) = F\(x). As a result we get Eqs. (19) and (20). 

Let us use the EM approximation. We neglect fluctu- 
ations around a mean value T(k, A) and use an approxi- 
mation 



F x (k,x) w e 



-ixT(k,X) 



(42) 



Then we get a self-consistent equation for the complex 
function T(k, A): 



T(fc,A) = £ 



P(ki\k) 



ki 



X + ie-Qa- l)T(fci,A) 



At A ^> 1 this equation has a solution 

T(k,X) = A -1 — iirXP (k\\ k) . 



(43) 



(44) 



This solution gives p(X) = 2k\P{k\) / |A| where k\ = 
A 2 + O(l) as before. It agrees with the result presented 
in Eq. (31) for uncorrelated graphs. One concludes that 
the short-range correlations between degrees of neighbor- 
ing vertices in the scale- free networks does not change the 
eigenvalue exponent 8 = 27 — 1. 



VII. SPECTRUM OF A TRANSITION MATRIX 

Let us consider random walks on a graph with the tran- 
sition probability l/k v of moving from a vertex v to any 
one of its neighbors. The transition matrix P then satis- 
fies 

Pw,v — Q>w,v /kv . (45) 
Clearly, for each vertex v 

J2 p W ,v = l- (46) 

w 

P is related with the Laplacian of the graph 

j 1 if v = u 

L„.,„ = < , (47) 



—<iv,w/\/k v ~k u ] otherwise 
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as follows 

P = D 1 / 2 (l-L)D~ 1 / 2 , (48) 

where D v w — S v _ w /k v . Therefore, if we know the den- 
sity p(X) of eigenvalues of P, we can find the density of 
eigenvalues of the Laplacian: pl(X) = p(l — A). 

We denote the eigenvalues of the matrix P by Ai ^ 
A2 ^ ... ^ Xn- The eigenfunction f v = k v corresponds 
to the largest eigenvalue Ai = 1. 

In order to calculate the spectrum of P we use the same 
method of random walks described in the Section II. The 
probability of one step is given by Eq. (45). We define 
the generating function Q v (z) and T^J (z^ 1 ) = TwJ{z) 
and obtain an exact equation which is similar to Eq. (13): 

1 



WV \ ) 



jZ ^ g Tg W (Z) 



(49) 



where g ~ w but g 7^ v. At z = X + ie, we get exact equa- 
tions for the function F\(x) = ^exp(— ixTw}(X + is)) 
and the density of the eigenvalues p(X): 



Fx(x) 



Jo yV 

(50) 

p(A) = -Rc^P(fc)fc / dye^F^y). (51) 

The function F\(x) = e~ lx is an exact solution of Eq. 
(50). This solution corresponds to the eigenvalue Ai = 1 
and gives the delta-peak 6(X — 1) in the density p(X). 
The second largest eigenvalue A2 is related to several im- 
portant graph invariants such as the diameter D of the 
graph, sec, for example, [9]: 



D(G) < 



log(7V - 1) 



(52) 



log(l/A 2 ) 

Here the diameter of a graph is the maximum distance 
between any two vertices of a given graph. 

In order to find the spectrum at A ^ A 2 we use the 
EM approach. We assume F\(x) w e - lxT W and get an 
equation for a complex function T(A) : 



(k) 

p(X) is given by 



kX + ie- (k- l)T(A) 



1 T 1 

T Im A-T(Ay 



(54) 



For completeness, we present the spectrum of the transi- 
tion matrix of a /c-regular tree: 



k sjA{k-l)/k 2 - X 2 
27 1 - A 2 



(55) 



which easily follows from Eqs. (53) a nd (54). The second 
eigenvalue is equal to A 2 = 2\Jk — 1/k. 



VIII. ANALYSIS OF SPECTRA 



Let us compare available spectra of classical random 
graphs and scale-free networks [13,14], empirical spectra 
of the Internet [10,11,16], and spectra of random tree-like 
graphs. 

At first we discuss spectra of adjacency matrices. The 
spectra were calculated in the framework of the EM ap- 
proach from Eqs. (23) and (24) for different degree dis- 
tributions P(k). Our results are represented in Figs. 1 
and 2. 

Classical random graphs. Classical random graphs 
have the Poisson degree distribution. The density of 
eigenvalues of the associated adjacency matrix has been 
obtained numerically in [13]. In Fig. 1 we display results 
of the numerical calculations and our results obtained 
within the EM approach. We found a good agreement 
in the whole range of eigenvalues. There are only some 
small differences in the region of small eigenvalues which 
may be explained by an inaccuracy of the EM approach 
in this range. In this region, the density p(X) has an ele- 
vated central part that differs noticeably from the semi- 
circular distribution. The spectrum also has a tiny tail 
given by Eq. (31) which can hardly be seen in Fig. 1, 
see for detail Section V and Refs. [21,22] 




X /<k> 



FIG. 1. Density of eigenvalues of the adjacency matri- 
ces of two networks, (i) The classical random graph (the 
Erdos-Renyi model) with the average degree (fc) = 10: the 
effective medium (EM) approach (the solid line) and numer- 
ical calculations for the graphs of 20 000 vertices [13] (the 
open circles), (ii) The scale- free random tree-like graph with 
7 = 3 and the smallest degree fco = 5: the EM approach 
(the dashed line), the improved EM approach, see the text 
(the dashed-dotted line). The results of the simulations of 
the Barabasi- Albert model of 7000 vertices [13] (the open 
squares). The semicircular law is shown by the thin solid 
line. 
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FIG. 2. Evolution of spectra of a random tree-like graph 
with the scale- free degree distribution for 7 = 5 and the small- 
est degree ko = 1, 2, 3 and 5. The panel a shows the spectra of 
the graphs with fco = 2, 3 and 5. The dotted line corresponds 
to the density of eigenvalues of an infinite chain. The dashed 
and thin solid lines present the spectrum of the k = 3 and 6 
regular Bethe lattices. The panel b shows the spectrum of a 
random uncorrelated graph having dead-end vertices with the 
probability P(l) = 0.3. The insert shows the behavior of the 
density of eigenvalues p(A) near the dip at P(l) = 0.1,0.2, 
and 0.3. The central peak is produced by localized states. 

Scale-free networks. Spectra of scale-free graphs with 
the degree distribution P(k) = P^k^ 1 differ strongly 
from the semicircular law [13,14]. The Barabasi- Albert 
model has a tree-like structure, the exponent 7 = 3 of 
the degree distribution, and negligibly weak correlations 
between degrees of the nearest neighbors [3]. Therefore, 
one can assume that the spectrum of a random tree-like 
graph can mimic well the spectrum of the model. In 
Fig. 1 we compare the spectrum of the random tree-like 
graph with 7 = 3 and the spectrum of the Barabasi- 
Albert model obtained from simulations [13]. The den- 
sity of states p(X) has a triangular-like form and demon- 
strates a power-law tail. There is only a noticeable de- 
viation of the EM results from the results of simulations 
[13] at small eigenvalues A. In order to improve the 
EM results, we used, as an ansatz, the distribution func- 



tion F\(x) = [1 + a(X)x 2 ]e~ lxT ^ instead of the function 
F\(x) = e~ lxT ( x \ In this case, there are two unknown 
complex functions a(A) and T(A) which were determined 
self-consistently from Eq. (19). 

Power-law tail. The power-law behavior of the den- 
sity of eigenvalues p{\) oc \~ s is an important feature of 
the spectrum of scale-free networks. The simulations [13] 
of the Barabasi- Albert model having the degree exponent 

7 = 3 revealed a power-law tail of the spectrum, with the 
eigenvalue exponent 5 w 5. Our prediction S = 27— 1 = 5 
is in agreement with the result of these simulations. 

The study of the topology of the Internet at the Au- 
tonomous System (AS) level revealed a power-law be- 
havior of eigenvalues of the associated adjacency ma- 
trix [10,11]. The degree distribution of the network 
has the exponent 7 sa 2.1 [11]. The eigenvalues \ of 
the Internet graph are proportional to the power of the 
rank i of an eigenvalues (starting with the largest eigen- 
value): A, oc i e with some exponent e. This leads to 
p(A) oc J2i S(X - Xi) oc A~ 1+1 / e . The Multi dataset ana- 
lyzed in [11] gave e w —0.447 and, hence, the eigenvalue 
exponent 1 — 1/e w 3.2. The Oregon dataset [11] gave 
e w -0.477, 1 - l/e w 3.1. 

Our results with 7 = 2.1 substituted, give the eigen- 
value exponent 5 = 27 — 1 3.2 in agreement with the 
results obtained from empirical data for this network. 
There are the following reasons for the agreement be- 
tween the theory for tree-like graphs and the data for 
the Internet. At first, although the average clustering 
coefficient of the Internet at AS level is about 0.2, the 
local clustering coefficient rapidly decreases with increas- 
ing the degree of a vertex [36]. In other words, the clos- 
est neighborhood of vertices with large numbers of con- 
nections is "tree-like". Recall that vertices with large 
numbers of connections determine the large-eigenvalue 
asymptotics of the spectrum. So, we believe that our 
results for the asymptotics of the spectra of tree- like net- 
works is also valid for the Internet and other networks 
with similar structure of connections. Secondly, the In- 
ternet is characterized by strong correlations between de- 
grees of neighboring vertices [33]. However, as we have 
shown in the Section VI, such short-range degree correla- 
tions do not affect the power-law behavior of eigenvalues. 

The study of the Internet topology [11] also revealed 
a correspondence between the large eigenvalues Xi and 
the degree kf. ki — Af. This result is in agreement with 
our theoretical prediction that it is the highly connected 
vertices with a degree about k\ w A 2 that produce the 
power-law tail p(A) oc A~ 5 . 

The calculations of the eigenvalues spectrum of the 
adjacency matrix of a pscudofractal graph with 7 = 
2.585 . . . [37] have revealed a power-law behavior with 

8 w 4.6. The effective medium approximation gives lower 
value 5 = 2-y—l w 4.2. The origin of the difference is not 
clear. One should note that the pseudofractal is a deter- 
ministically growing graph with a very large clustering 
coefficient C = 4/5 and, what is especially important, 
with long-range correlations between degrees of vertices. 
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FIG. 3. Local configurations which produce localized states 
with zero-eigenvalue (the central peak) in the spectrum of a 
random graph (see the text). 

Weakly connected nodes. Let us study the influence 
of weakly connected vertices with degrees 1 $5 fc 5 on 
the spectra of random tree-like graphs with the degree 
distribution P(k) = P^k^ 1 . In Fig. 2a and 2b we rep- 
resent the evolution of the spectrum of the network with 
7 = 5, when the smallest degree fc decreases from 5 to 
1. The spectra were calculated in the framework of the 
EM approximation. Similar results are obtained at dif- 
ferent 7. For ko ^ 4, two peaks at non-zero eigenvalues 
emerge in the density of states p(A). In order to under- 
stand an origin of the peaks one can note that for this 
degree distribution the average degree (k) is close to ko. 
For example, at fco = 3 we have (k) = 3.49. Therefore, 
in this network, the probability to find a vertex having 
three links is larger than the probability to find a vertex 
with a degree k ^ 4. There are large parts of the network 
which have a local k — 3-regular structure. In Fig. 2a 
we show a density of eigenvalues of an infinite k = 3 
regular Bethe lattice [see Eq. (16) at k = 3]. At small 
eigenvalues, the density of the regular tree fits well the 
density of the random network. At large A, the density 
of eigenvalues demonstrates a power-law behavior with 
the exponent S = 27 — 1. 

In the case fco = 2 we have (fc) = 2.23. This network 
contains long chains which connect vertices with degrees 
fc ^ 3. In Fig. 2a we display the density of eigenvalues 
of an infinite chain (sec Eq. (16) at fc = 2). At small 
eigenvalues this density of eigenvalues fits well the den- 
sity of eigenvalues of the random network. Therefore, 
it is the vertices with small degrees that are responsible 
for the formation of density p(A) of networks at small 
eigenvalues. 

Dead-end vertices. Let us investigate the effect of dead- 
end vertices on the spectra of random tree-like graphs 
with different degree distributions. Fig. 2b shows a spec- 
trum of a scale- free network with 7 = 5 and the probabil- 
ity of dead-end vertices P(l) = 0.3. The EM approxima- 
tion is used. The spectrum has a flat part and two peaks 
at moderate eigenvalues. As we have shown above, this 
(intermediate) part of the spectrum is formed mainly by 
the vertices with degree fc = 2 and 3. The emergency of 
a dip at zero is a new feature of the spectrum. In fact, 



there is a gap in the spectrum obtained in the in the 
framework of the EM approach. The width of the gap 
increases with increasing P(l). One can see this in the 
insert on the Fig. 2b. The dead-end vertices also produce 
a delta peak at A = 0. The central peak corresponds to 
localized eigenstates. 



Note that the appearance of the central peak and a dip 
is a general phenomena in random networks with dead- 
end vertices. We also observed this effect in the classical 
random graphs. Spectral analysis of the Internet topol- 
ogy on the AS level revealed a central peak with a high 
multiplicity [16]. Thus the conjecture that localized and 
extended states are separated in energy may well hold 
in complex networks. A similar spectra was observed 
in many random systems, for example, in a binary alloy 
[38]. In order to estimate the height of the delta peak 
it is necessary to take into account all localized states. 
Unfortunately, so far this is an unsolved analytical prob- 
lem [16]. In Fig. 3 we show local parts of a network, 
which produce localized states. One can prove that con- 
figurations with two and more dead-end vertices, see Fig. 
3a, produce eigenstates with A = . The corresponding 
eigenvectors have non-zero components only at the dead- 
end vertices [16,39]. Fig. 36 shows another configuration 
which produces an eigenstate with the eigenvalue A = 0. 
A corresponding eigenvector is localized at vertices 0, 1 
and 2. 



Finite-size effects. In the present paper we studied the 
spectral properties of infinite random tree-like graphs. 
Numerical studies of large but finite random trees demon- 
strate that the spectrum of a finite tree consists, speak- 
ing in general terms, of a continuous component and an 
infinity of delta peaks. The components correspond to 
extended and localized states, respectively [17]. There is 
a hole around each delta peak in the spectrum. A finite 
regular tree has a spectral distribution function which 
looks like a singular Cantor function [39]. These results 
demonstrate that finite size effects in spectra may be very 
strong. In particular, the finite size of a network deter- 
mines the largest eigenvalue in its spectrum. As was 
estimated in the Section V, the largest eigenvalue of the 
adjacency matrix associated with a scale-free graph is of 
the order of fc^ 2 = fcfjV 1 / 2 ^- 1 ). 



Spectrum of the transition matrix. In Fig. 4 we repre- 
sent a spectrum of the transition matrix P defined by Eq. 
(45) for a tree-like graph with the scale-free degree distri- 
bution P(k) cc fc~ 7 at large degrees fc > 5. The spectrum 
was calculated from Eqs. (53) and (54) with the degree 
exponent 7 = 2.1 and the probabilities P(l), P(2), P(3) 
and P(4) taken from empirical degree distribution of the 
Internet at the AS level [36]. 
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IX. CONCLUSIONS 




X 

FIG. 4. Density of eigenvalues of the transition matrix P 
defined by Eq. (45). (i) The Internet data (the solid cir- 
cles) and the results of simulations of a random scale-free 
networks (the open squares) from Ref. [15]. (ii) The spec- 
trum of the Internet extracted from Ref. [16] (the stars) . (iii) 
Our calculations (the solid line) with the degree distribution 
P(k) = Ak- X1 for k > 4, P(l) = 0.358, P(2) = 0.4 
P(3) = 0.12. These parameters are taken from Ref. [36]. The 
height of the central peak was estimated from Ref. [16]. The 
insert shows the spectra of the transition matrix of a random 
tree-like graph with excluded dead-end vertices: (i) a random 
tree-like graph with a scale-free degree distribution, 7 = 3, 
the smallest degree feo = 5, and (k) = 9.06 (the dashed line), 
(ii) A classical random graph with (k) = 9.06 (the dotted 
line), (iii) A k— regular Bethe lattice with k — 9 (the solid 
line) . 



The spectrum lies in the range |A| ^ A2 < 1. In Fig. 
4 we compare our results with the spectrum of the tran- 
sition matrix P of the Internet obtained in [15,16]. Un- 
fortunately, the data [15,16] are too scattered to make a 
detailed comparison with our results. Nevertheless, one 
can see that the spectrum of P of the tree-like graph re- 
produces satisfactory the general peculiarities of the real 
spectrum. Namely, the spectra have a wide dip at zero 
eigenvalue and a central delta-peak [16]. The multiplicity 
of the zero eigenvalue have been estimated in [16]. For a 
detailed comparison between the spectra, correlations in 
the Internet must also be taken into account. 

In order to reveal &n effect of dead-end vertices we cal- 
culated spectra of P on a random tree-like graph with 
the Poisson and the scale-free degree distributions 7 = 3 
in the case when dead-end vertices are excluded, that is 
P(l) = 0, and (k) = 9.06. These spectra are displayed 
in the insert on the Fig. 4. In the whole range of eigen- 
values these spectra are very close to the spectrum of a 
fc— regular Bethe lattice with the degree k — 9. These cal- 
culations confirm the fact that it is the dead-end vertices 
that produce the dip in the spectrum of the Internet. 



In this paper we have studied spectra of the adjacency 
and transition matrices of random uncorrelated and cor- 
related tree-like complex networks. We have derived ex- 
act equations which describe the spectrum of random 
tree- like graphs, and proposed a simple approximate solu- 
tion in the framework of the effective medium approach. 
Our study confirms that spectra of scale-free networks 
as well as the spectra of classical random graphs do not 
satisfy the Wigner law. 

We have demonstrated that the appearance of a tail of 
the density of the eigenvalues of sparse random matrices 
is a general phenomenon. The spectra of classical random 
graphs (the Erdos-Renyi model) have a rapidly decreas- 
ing tail. Scale-free networks demonstrate a power-law 
behavior of the density of eigenvalues p(A) oc |A| S . We 
have found a simple relationship between the degree ex- 
ponent 7 and the eigenvalue exponent 5: 5 = 27 — 1 . We 
have shown that correlations between degrees of neigh- 
boring vertices do not affect the power-law behavior of 
eigenvalues. Comparison with the available results of 
the simulations of the Barabasi-Albert model and the 
analysis of the Internet at the Autonomous System level 
shows that this relationship is valid for these networks. 
We found that large eigenvalues A > 1 are produced by 
highly connected vertices with a degree k w A 2 . 

Many real-world scale-free networks demonstrate 
short-range correlations between vertices [35,40] and a 
decrease of a local clustering coefficient with increas- 
ing degree of a vertex. Therefore, the relationship S = 
27 — 1 between the degree-distribution exponent 7 and 
the eigenvalue exponent S may also be valid for these 
networks. We can conclude that the power-law behav- 
ior p(A) oc \~ s is a general property of real scale-free 
networks. 

Weakly connected vertices form the spectrum at small 
eigenvalues. Dead-end vertices play a very special role. 
They produce localized eigenstates with A = (the cen- 
tral peak). They also produce a dip in the spectrum 
around the central peak. In conclusion, we believe that 
our general results for the spectra of tree-like random 
graphs are also valid for many real- world networks with 
a tree-like local structure and short-range degree corre- 
lations. 
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Note added. — After we have finished our work we have 
learned about a recent mathematical paper Ref. [41], 
where large eigenvalues of spectra of complex random 
graphs were calculated. The statistical ensemble of 
graphs, which was considered in that paper, essentially 
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differs from that of our paper and has a different cutoff 
of the degree distribution, but the asymptotics of spectra 
agree in many cases. 
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